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Abstract 

■ We study the LHC discovery potential for heavy Major ana neutrino singlets 
O ' in the process pp — > — > £~^N i'^i'^jj = e, /i) plus its charge conjugate. 

' With a fast detector simulation we show that backgrounds involving two like-sign 

' charged leptons are not negligible and, moreover, they cannot be eliminated with 

I simple sequential kinematical cuts. Using a likelihood analysis it is shown that, 

I for heavy neutrinos coupling only to the muon, LHC has 5a sensitivity for masses 

^ I up to 200 GeV in the final state n^ji^jj. This reduction in sensitivity, compared 

J> \ to previous parton-level estimates, is driven by the ~ 10^ — 10^ times larger 

^ \ background. Limits are also provided for e^e^jj and e^^^jj final states, as 

I well as for Tevatron. For heavy Dirac neutrinos the prospects are worse because 

Q I backgrounds involving two opposite charge leptons are much larger. For this case, 

I we study the observability of the lepton flavour violating signal e^/i^jj. As a 

I by-product of our analysis, heavy neutrino production has been implemented 

Qh! within the ALPGEN framework. 

^ ■ 

> ; 1 Introduction 

X: 

■ Large hadron colliders involve strong interacting particles as initial states, giving rise 
to huge hadronic cross sections. The large luminosities expected will also provide quite 
large electroweak signals, with for instance 1.6 x 10^° (4 x 10^) W bosons at LHC 
(Tevatron) for a luminosity of 100 (2) fb~^. Therefore, these colliders can be used 
for precise studies of the leptonic sector as well, and in particular they can produce 
new heavy neutrinos at an observable level, or improve present limits on their masses 
and mixings [1-4] (see Ref. [5] for a review). These new fermions transform trivially 
under the gauge symmetry group of the Standard Model (SM), and in the absence of 
other interactions they are produced and decay only through their mixing with the 
SM leptons. With new interactions, like for instance in left-right models [6], heavy 
neutrinos can be produced by gauge couplings unsuppressed by small mixing angles. 
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yielding larger cross sections and implying a much higher collider discovery reach [7-10]. 
Heavy neutrinos could also be copiously produced in pairs through the exchange of a 
relatively light Z' boson [11]. In these scenarios, however, the observation of the new 
interactions could be more interesting than the existence of new heavy neutrinos. 

We will concentrate on the first possibility and neglect other new production mechan- 
ims, taking a conservative approach. In this case, for example, it has been claimed by 
looking at the lepton number violating (LNV) AL = 2 process — > /i^yU^jj that LHC 
will be sensitive to heavy Majorana neutrinos with masses mjy up to 400 GeV, whereas 
Tevatron is sensitive to masses up to 150 GeV [2,4]. However, as we will show, taking 
into account the actual backgrounds these limits are far from being realistic. In par- 
ticular, backgrounds involving h quarks, including for instance tinj (with nj standing 
for n = 0, 1, 2, . . . additional jets), are two orders of magnitude larger than previously 
estimated. Moreover, in the region m^v < M^i/ the largest and irreducible background 
is 66nj, by far dominant but overlooked in previous parton-level analyses [4]. In this 
work we make a detailed study, at the level of fast simulation, of the LHC sensitivity 
to Majorana neutrinos in the process /i^/i^jj, which is the cleanest final state, 

for both rriN > My/ and mjq < M^r. We also study the processes pp — > e^e^jj and 
pp —>■ e^^'^jj for which the sensitivity is slightly worse. Heavy Dirac neutrinos do not 
produce LNV signals and then their observation is much more difficult. As an example, 
we examine the lepton flavour violating (LFV) signal e^/i^jj, produced by a heavy 
Dirac neutrino coupling to the electron and muon. 

The generation of heavy neutrino signals has been implemented in the ALPGEN [12] 
framework, including the process studied here as well as other final states. In the 
following, after making precise our assumptions and notation in section [21 we describe 
the implementation of heavy neutrino production in ALPGEN in sectional We present 
our detailed results in section HI where we will eventually find that heavy neutrinos can 
be discovered up to masses of the order of 200 GeV, and that for N lighter than the 
W boson its mixing can be probed at the 10~^ level (for a "reference" mass mjv = 60 
GeV). These figures are much less optimistic than in previous literature. Estimates for 
Tevatron are given in section [5], and our conclusions are drawn in section [6l In two 
appendices we detail the evaluation of the bbrij background and the heavy neutrino 
mass reconstruction, respectively. 
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2 Heavy neutrino interactions 



Our assumptions and notation are reviewed in more detail in Ref. [5] (see also Refs. [13, 
14]). The SM is only extended with heavy neutrino singlets Nj, which are assumed 
to have masses of the order of the electroweak scale, up to few hundreds of GeV. We 
concentrate on the lightest one, assuming for simplicity that the other extra neutrinos 
are heavy enough to neglect possible interference effects. The new heavy neutrino 
(where we suppress the unnecessary subindex) can have Dirac character, what requires 
the addition of at least two singlets, or Majorana, in which case {N^Y = CNj^ = Nji 
and lepton number is violated. In either case it is produced and decays through its 
mixing with the light leptons, which is described by the interaction Lagrangian (in 
standard notation) 



Cw = -^{iYVeNPLNW^ + NYV;j,PLiWl) 
9 



Cz = -^{ua^V,MPLN + Nrv;^PLy,)z,, 

gniN 



'^H = -^!-^ {viV,NPRN + NV;^PLV,)H. (1) 

The SM Lagrangian remains unchanged in the limit of small mixing angles Vin-, ^ = 
e, /i, r (which is the actual case), up to very small corrections 0{V^). Neutral couplings 
involving two heavy neutrinos are also of order V^. The heavy neutrino mass rriN- joins 
two different bispinors in the Dirac case and the same one in the Majorana case. Heavy 
neutrino decays are given by their interactions in Eqs. ([T]): N ^ W~^i~, N Zu, 
N Hu, plus — i> W~i~^ for a heavy Majorana neutrino. For rriN < My/ all these 
decays produce three body final states, mediated by off-shell W , Z or H bosons. The 
total width for a Majorana neutrino is twice larger than for a Dirac one with the same 
couplings [5,15-17]. 

As it is apparent from Eqs. ([T]), heavy neutrino signals are proportional to the 
neutrino mixing with the SM leptons Vin- Limits on these matrix elements have been 
extensively discussed in previous literature, and we quote here only the main results. 
Low-energy data constrain the quantities 

3 n 

= See' - ^ Viuye*u, = ^ ^iN^Vf^N, ■ (2) 
i=i j=i 

A global fit to tree level processes involving light neutrinos as external states gives 
[18,19], 

^ee < 0.0054 , < 0.0096 , Q^r < 0.016 (3) 
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at 90% confidence level (CL). Note that a global fit without the unitarity bounds 
implies Qee < 0.012 [18]. Additionally, for Majorana neutrinos coupling to the electron 
the experimental bound on neutrinoless double beta decay requires [20] 



2 



< 5 X 10"* GeV"^ 



(4) 



If V'eTVj saturate flee in Eq. ([3]), this limit can be satisfied either demanding that rriNj 
are large enough, beyond the TeV scale [21] and then beyond LHC reach, or that 
there is a cancellation among the different terms in Eq. (jlj), as may happen in definite 
models [22], in particular for (quasi)Dirac neutrinos. 

Flavour changing neutral processes further restrict Qw. The new contributions, and 
then the bounds, depend on the heavy neutrino masses. In the limit m^, ^ ^ 







they imply [24] 



l^e^l < 0.0001 



\Q,J < 0.01 



\%t\ < 0.01 



(5) 



Except in the case of flefi, for which experimental constraints on lepton flavour violation 
are rather stringent, these limits are similar to the limits on the diagonal elements. An 
important difference, however, is that (partial) cancellations among loop contributions 
of different heavy neutrinos may be at work [25]. Cancellations with other new physics 
contributions are also possible. Since we are interested in determining the heavy neu- 
trino discovery potential and the direct limits on neutrino masses and mixings which 
can be eventually established, we must consider the largest possible neutrino mixings, 
although they may require model dependent cancellations or fine-tuning. 



3 Heavy neutrino production with ALPGEN 

For the signal event generation we have extended ALPGEN [12] with heavy neutrino 
production. This Monte Carlo generator evaluates tree level SM processes and provides 
unweighted events suitable for simulation. A simple way of including new processes 
taking advantage of the ALPGEN framework is to provide the corresponding squared 
amplitudes decomposed as a sum over the different colour structures. In the case of 
heavy neutrinos this is trivial because there is only one term. This method requires to 
evaluate from the beginning the squared amplitudes for the processes one is interested 

^When ViNj > Mw/rriN^ the non-decoupling terms in the amphtude, proportional to 
'^iNj^Nj/^w^ cannot be neglected because they dominate over the V^ff. terms [23]. 
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in, what is done using HELAS [26]. An alternative possibility which gives more flexi- 
bility for future applications is to implement the new vertices at the same level as the 
SM ones, what is quite more involved. 

We have restricted ourselves to single heavy neutrino production. Pair production 
is suppressed by an extra V"^ mixing factor and by the larger center of mass energy 
required, what implies smaller PDFs and more suppressed s-channel propagators. Sin- 
gle heavy neutrino production can proceed through s-channel W, Z or H exchange. 
The flrst two production mechanisms have been implemented in ALPGEN for the var- 
ious possible flnal states given by the heavy neutrino decays W^i^, N —>■ Zui, 
N Hui with £ = e, /i, r, and for both Dirac or Majorana A^. In the case rriN < My/ 
all decays are three-body, and mediated by off-shell W ^ Z or H. The transition from 
two-body to three-body decays on the M^y, Mz and thresholds is smooth, since 
the calculation of matrix elements and the A^ width are done for off-shell intermediate 
bosons. Two approximations are made, however. The small mixing of heavy neu- 
trinos with charged leptons implies that their production is dominated by diagrams 
with A^ on-shell, like those shown in Fig. [H with a pole enhancement factor, and that 
non-resonant diagrams are negligible. (Additionally, to isolate heavy neutrino signals 
from the background one expects that the heavy neutrino mass will have to be recon- 
structed to some extent.) Then, the only diagrams included are the resonant ones. In 
the calculation we also neglect hght fermion masses except for the bottom quark. 



Figure 1: Feynman diagrams for the process qq' i^N, followed by LNV decay 
A^ — > i^'^^'^W" (a) and lepton number conserving (LNC) decay A^ — > i^'^^'W^ (b). The 
diagrams for the charge conjugate processes are similar. 

Generator- level results are presented in Fig. [2] for LHC and Tevatron in the relevant 
mass ranges. Sohd lines correspond to the total fiN cross sections for |V^Ar| = 0.098, 
VeN = VrN = 0. The dashed lines are the cross sections for the flnal state /i^/x^jj. 




(a) 



(b) 
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which is the cleanest one. The dotted hnes are the same but with kinematical cuts 
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> 0.4, 
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< 2.5, 
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> 10 GeV, 


\V''\ 


< 2.5, 







(6) 



included to reproduce roughly the acceptance of the detector and give approximately 
the "effective" size of the observable signal. Of course, the correct procedure is to per- 
form a simulation, as we do in next section, but for illustrative purposes we include the 
cross-sections after cuts. In particular, they clearly show that although for ttln < Mw 
the total cross sections grow several orders of magnitude, both at LHC and Tevatron, 
partons tend to be produced with low transverse momenta (the two muons and two 
quarks result from the decay of an on-shell W), making the observable signal much 
smaller. These results are in agreement with those previously obtained in Ref. [4]. 





Tevatron 



' 40 50 60 70 80 90 100 110 120 130 140 



Figure 2: Cross sections for heavy neutrino production at LHC (left) and Tevatron 
(right), as a function of the heavy neutrino mass, for |V"^Ar| = 0.098. The solid lines 
correspond to total [xN cross section, the dashed lines include the decay to like-sign 
muons and the dotted lines are the same but including the kinematical cuts in Eq. 



4 Di-lepton signals at LHC 

The most interesting scenario for LHC is when the heavy neutrino has Majorana nature 
and couples only to the muon, so that it produces a final state /i^/i^jj with two 
same sign muons and at least two jets. Since this LNV signal has sometimes been 
considered [2,4] to be almost background free (more realistic background estimates are 
given in Ref. [27]), a detailed discussion of the actual backgrounds is worthwhile. A 
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first group of processes involves tlie production of additional leptons, either neutrinos or 
charged leptons (which may be missed in the detector). The main ones are W^W^nj 
and W^Znj. We point out that not only the processes with n = 2 contribute: processes 
with n < 2 are backgrounds due to the appearance of extra jets from pile-up, and 
processes with n > 2 cannot be cleanly removed because of pile-up on the signal. A 
second group includes final states with b and/or b quarks, like ttnj, with semileptonic 
decay of the ti pair, and Wbbnj, with W decaying leptonically. In these cases the 
additional like-sign muon results from the decay of a 6 or 6 quark. Only a tiny fraction 
of such decays produce isolated muons with sufficiently high transverse momentum. 
But, since the tinj and Wbbnj cross sections are so large, these backgrounds are also 
much larger than backgrounds with two weak gauge bosons. Finally, bbnj production is 
several orders of magnitude larger than all processes mentioned above, but the produced 
muons have small pt and invariant mass in this case. Then, in general it might be 
eliminated with suitable high-px' cuts on charged leptons [28] (see section I^TTl) . but for 
rriN < M\Y the heavy neutrino signal is also characterised by very small transverse 
momenta (see section and this background turns out to be the dominant one. 
The same applies for ccnj, but with the difference that c quark decays produce isolated 
charged leptons much less often than b decays. 

Other LNV signals produced by heavy neutrinos are e^e^jj and e^fi^jj. They 
have the same SM backgrounds but with one important difference: b decays produce 
"apparently isolated" electrons more often than muons, because electrons are detected 
in the calorimeter while muons travel to the muon chamber. Hence, the correspond- 
ing backgrounds tinj, bbnj e^e^X/e^/i^X are larger than the ones involving only 
muons. A precise evaluation of these backgrounds, optimising the criteria for electron 
isolation, seems to require a full simulation of the detector. The limits provided in these 
cases must be regarded with some caution in this respect, and should be confirmed with 
a full detector simulation. 

We have generated the signal and backgrounds using ALPGEN and passing them 
through PYTHIA 6.4 [29] with the MLM prescription [30] to avoid double counting of 
jet radiation. A fast simulation of the ATLAS detector [31] has been performed. For 
the signal and all backgrounds except bbnj and ccnj the number of simulated events 
corresponds to at least 10 times the luminosity considered (which is 30 fb"'^), so as to 
reduce statistical fluctuations, and the number of events is scaled accordingly. For bbnj 
and ccnj the luminosity simulated is 0.075 fb^^. Their evaluation is further discussed 
in appendix [Al It must also be noted that in the signal simulation all W decays in 
pp iN WW are included. Leptonic W decays give an extra ~ 20% contribution 
to di-lepton final states when the charged lepton from the W decay is missed, or when 
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W decays to rz/ and the tau lepton decays hadronically. 
4.1 (-^(-^33 production for rriN > Mw 

In this mass region we take the reference values rriN = 150 GeV and (a) V^n = 0.098, 
VeN = Ktv = 0; (b) Ktv = 0.073, V^^ = Ktv = 0; (c) Ktv = 0.073, V^^ = 0.098, 
VrN = 0. The pre-selection criteria used for our analysis are: 

(i) two like-sign isolated charged leptons with pseudorapidity \ri^\ < 2.5 and trans- 
verse momentum larger than 10 GeV (muons) or 15 GeV (electrons), and no 
additional isolated charged leptons; 

(ii) no additional non-isolated muons; 

(iii) two jets with \r]^ < 2.5 and > 20 GeV. 

We point out that for fi'^fi'^jj final states the requirement (ii) reduces the backgrounds 
involving Z bosons by almost a factor of two, and thus proves to be quite useful. 
The number of events at LHC for 30 fb~^ after pre-selection cuts is given in Table [H 
Additional backgrounds such as tb, titi, tibb, Ztinj, WWZnj, WZZnj and ZZZnj 
are smaller and we do not show them, but they are included in the estimation of the 
signal significance below. The number of like-sign dimuon events from ccnj displayed 
between parentheses corresponds to an estimation, because no /i^/i^X events are found 
in the sample simulated (more details can be found in appendix |A]). We also note that 
the higher threshold for electrons contributes to the difference between the numbers 
of e^e^jj and /x^e^jj events, which are expected to be similar in some cases, for 
example for tinj. 

Let us concentrate on fi'^fi'^jj final states. The fast simulation shows that SM 
backgrounds are about two orders of magnitude larger than previously estimated (three 
orders if we include bbnj). Moreover, they cannot be sufficiently suppressed with 
respect to the heavy neutrino signal using simple cuts. Some obvious discriminating 
variables are: 

• The missing momentum /j-. It is smaller for the signal because it does not have 
neutrinos in the final state, but nonzero due to energy mismeasurement in the 
detector. 

• The separation between the muon with smallest Pt (we label the two muons as 
/ii, fi2, by decreasing transverse momentum) and the closest jet, AR^^j. For 
backgrounds involving high-pr b quarks this separation tends to be rather small. 
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Table 1: Number of i^i^jj events at LHC for 30 ih-\ at the pre-selection and selection 
levels. The heavy neutrino signal is evaluated assuming m^r = 150 GeV and coupling 
(a) to the muon, V^n = 0.098; (b) to the electron, VeN = 0.073; (c) to both, VeN = 0.073 
and Vf,N = 0.098. 

• The transverse momentum of the two muons, pt^^ and p^^, respectively. In par- 
ticular is a good discriminant against backgrounds from b quarks, because 
these typically have one muon with small px- 

These variables are plotted in Fig.[3]for the /i^/U^jj signal and the backgrounds grouped 
in three classes with common features: (a) bbnj, where both muons come from b quark 
decays (the contribution of ccnj is negligible); (b) tinj, tj and W/ Zbbnj, where one 
muon comes from a b quark; (c) backgrounds where both muons come from W/Z decays 
(mainly di-boson and tri-boson production). Kinematical cuts on the variables listed 
above do not render the /x^/i^ jj final state "background free" , as it is apparent from the 
plots (and we have explicitly checked). Indeed, for the large background cross sections 
in Table [H the overlapping regions contain a large number of background events, and 
they can be eliminated only by severely reducing the signal. However, a likelihood 
analysis using these and further variables can efficiently reduce the background. The 
additional variables are shown in Fig. |3J 
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Figure 3: Normalised distributions of several discriminating variables for the /^^/x^jj 
signal with mjv = 150 GeV and its backgrounds (see the text). 



• The invariant mass rrijj of the two jets with largest transverse momentum, which 
for the signal are assumed to originate from the W hadronic decay, and the 
invariant mass of fi2 (the muon with lowest pr) and these two jets, mw^i2- (Further 
details about the W and mass reconstruction can be found in appendix [Bl ) 
An important observation in this case is that in backgrounds involving b quarks 
this muon typically has a small px, displacing the background peaks to lower 
invariant masses. 

• The invariant mass of the two muons. 

• The separation between the muon with largest pr and the closest jet, AR^^^j. 

• The number of 6-tagged jets A^^ and jet multiplicities Nj. Especially the former 
helps to separate the backgrounds involving b quarks because they often have b- 
tagged jets. In this fast simulation analysis we have fixed the 6-tagging efficiency 
to 60%, but in a full simulation the b tag probability can be included in the 
likelihood function, improving the discriminating power of this variable. 
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Figure 4: Normalised distributions of several discriminating variables for the fi^^^jj 
signal with mjq = 150 GeV and its backgrounds (see the text). 
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• The transverse momenta of the two jets with largest pt, p™^^ and p^^^"^ respec- 
tively. 

These variables are not suited for performing kinematical cuts but greatly improve the 
discriminating power of a likelihood function. The resulting log-likelihood function is 
also shown in Fig. El where we distinguish four likelihood classes as in the previous 
figures: the signal, bbnj, backgrounds with one muon from h decays, and backgrounds 
with both muons from W/Z decays. 



Figure 5: Log-likelihood function for the ^i^^i^jj signal with mjv = 150 GeV and its 
backgrounds. 

The probability distributions built for fi^fi^jj final states are used for e'^e^jj and 
[M^e^ii as well. As selection criteria we require log^^ Ls / Lb > 1-4 for fi'^fi'^jj and 
logi^Ls/LB > 2.5 for e^e^jj and ii^e^jj final states, respectively, and that at least 
one of the two heavy neutrino mass assignments mvi/^i, '^nw^l2 is between 130 and 
170 Gevi The number of events surviving these cuts can be read on the right part 
of Table [H As it is apparent, the hkelihood analysis is quite effective in suppressing 
backgrounds, especially bbnj, tinj and W/Zbbnj. The resulting statistical significance 
for the heavy neutrino signals are collected in Table [21 assuming a "reference" 20% 
systematic uncertainty in the backgrounds (which still has to be precisely evaluated in 
a dedicated study). The limits on heavy neutrino masses and couplings depend on the 
light lepton they are coupled to. We can consider two extreme cases: 

(a) A 150 GeV heavy neutrino coupling only to the muon can be discovered for 

^The latter requirement assumes a previous knowledge of rriN- In the same way, the signal distri- 
butions for the likelihood analysis must be built for a fixed toa? value. Thus, experimental searches 
must be performed by comparing data with Monte-Carlo samples generated for different values of 
rriN. This procedure, although more involved than a search with generic cuts, provides much better 
sensitivity. 
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Table 2: Statistical significance of the heavy neutrino signals in the different channels, 
for a mass m^v = 150 GeV and coupling (a) to the muon, V^at = 0.098; (b) to the 
electron, VeN = 0.073; (c) to both, VeN = 0.073 and V^n = 0.098. 

mixings |V^Ar| > 0.054, and if no background excess is found the limits |V^Arp < 
0.97 (1.2) X 10^^ can be set at 90% (95%) CL, improving the ones from low 
energy processes (see section |2l) by a factor of 10. Heavy neutrino masses up to 
200 GeV can be observed with 5a at the LHC for V^^n = 0.098. 

(b) A 150 GeV heavy neutrino coupling only to the electron can be discovered for 
mixings \ VeN\ > 0.080 (excluded by the limits in section [2]), but if no background 
excess is found the limits |VeArp < 2.1 (2.5) x 10~^, which are slightly better 
than the one derived from Eq. ([3]), can be set at 90% (95%) CL. Heavy neutrino 
masses up to 145 GeV can be observed with 5a at the LHC for V^n = 0.073. 

For a heavy neutrino coupling to the electron and muon the limits depend on both 
couplings as well as on its mass. The combined limits for m^v = 150 GeV are displayed 
in Fig. O Except in the regions with VeN ~ or V^jy ~ 0, the indirect limit from /i — e 
LEV processes, also shown in this plot, is much more restrictive. 

These limits can be considered conservative in the sense that only the lowest-order 
signal contribution (without hard extra jets at the partonic level) has been included, 
and further signal contributions iNnj should improve the heavy neutrino observability. 
If the Higgs is heavier than 120 GeV the branching ratios Br(A^ — > W£) will increase 
as well. We also stress again that in the e^e^jj and fi^e^jj channels the evaluation of 
ttnj and other backgrounds with isolated electrons from b quarks must be confirmed 
with a full simulation, with an eventual optimisation of the isolation criteria. This is 
beyond the scope of the present work. 

It is worth explaining here in more detail why our results are much more pes- 
simistic than previous ones. With this purpose, we apply to signal and backgrounds 
the sequential kinematical cuts in Ref. [4]: 

• Missing energy < 25 GeV. 

• Lego-plot separation AR^j > 0.5. 
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Figure 6: Combined limits on VeN and V^at, for at = and rriN = 150 GeV. The red 
areas represent the 90% CL hmits if no signal is observed. The white areas correspond 
to the region where a combined statistical significance of 5cr or larger is achieved. The 
indirect limit from n — e LFV processes is also shown. 

• Dijet invariant mass 60 GeV < rrijj < 100 GeV, where the two jets are expected 
to come from the W boson in the case of the signal. 

The number of events for the signal and main backgrounds after these cuts are gathered 
in the left column of Table [3] (we do not show smaller backgrounds for brevity) . For 
tun = 150 GeV and V^^n = 0.098 the signal cross section is reduced to 1.7 fb, to be 
compared to ~ 2.2 fb in Ref. [4]. But our total background cross section after cuts 
amounts to 44 fb, while their estimate is of 0.04 fb. This difference by a factor of 1000 
arises mainly from the bbnj background, overlooked before, which is by far the largest 
one. But even if bbnj is not taken into account, the background cross section ~ 0.88 
fb is 20 times larger, due to: (i) ttnj, which was assumed negligible after cuts, and 
W/Zbbnj, also overlooked; (ii) the WZnj background, because parton- level analyses 
underestimate the probability of missing a lepton and thus its contribution; (iii) pile-up, 
which makes lower order processes {n < 2) contribute. All these backgrounds, collected 
in Table [3l can be compared to WWW, which was found to be the main background 
before. The resulting statistical significance of the signal, ignoring systematic errors, 
is S/^/B = 1.41cr for 30 fb~^, far from the ~ 30cr previously estimated. (If one makes 
the more realistic assumption that systematic errors are of order 20%, as we do in this 
work, then the statistical significance is further reduced to 0.19o".) It could be argued 
that the cuts in the previous list might be strengthened in order to further reduce the 
backgrounds. But this would be at the cost of reducing the signal as well. On the 
other hand, additional cuts on lepton transverse momenta can be introduced to reduce 
bbnj and tinj. Requiring that one charged lepton has px > 30 GeV and the other 
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one pt > 20 GeV, the signal is hardly affected while bbnj is essentially eliminated, as 
it is shown in the second column of Table El The statistical significance in this case 
is S/\fB = 14.1cr (ignoring systematic errors) or 12.1a (with 20% systematics) . We 
emphasise that, as it can be observed by comparing Tables [1] and [3], a probabilistic 
analysis is much more powerful in this case than a standard one based on cuts. But 
at any rate recovering parton-level estimates for the signal significance seems hardly 
possible. 

Sequential cuts I Sequential cuts II 



N (fi) 51.3 44.0 

bbnj 1293 2.7 

tinj 15.3 1.4 

Wbbnj 3.6 0.2 

Wttnj 0.7 0.7 

Zbbnj 0.9 0.0 

WWnj 0.5 0.5 

WZnj 4.1 2.9 

WWWnj 1.1 0.9 



Table 3: Number of fi^fJ^^jj events at LHC for 30 fb~^, after the kinematical cuts in 
Ref. [4] (first column) and with additional requirements (second column, see the text). 
The heavy neutrino signal is evaluated assuming rripq = 150 GeV and V^^ = 0.098. 

Finally, we would like to note that we have not addressed the observability of heavy 
neutrino signals in r lepton decay channels because they are expected to have much 
worse sensitivity. For hadronic r decays the charge of the decaying lepton seems rather 
difficult to determine, hence backgrounds from top pair and Z production will be huge 
(see also section 1^751 below). For leptonic decays r — > iv^vi, i = e.ii, not only the 
branching ratios are smaller, but also the signal has final state neutrinos and thus the 
discriminating power of against di-boson and tri-boson backgrounds is much worse. 

4.2 i'^i'^jj production for jun < Mw 

In this mass region we take the reference values mjv = 60 GeV and (a) V^n = 0.01, 
VeN = = 0; (b) Kiv = 0.01, V^N = Kiv = 0; (c) V^n = 0.01, V^n = 0.01, 
VrN = 0- The pre-selection criteria are the same as before. The likelihood analysis 
is performed distinguishing four classes: the signal, bbnj, backgrounds with one muon 
from b decays, and backgrounds with both muons from W/Z decays. The relevant 
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variables are depicted in Figs. [7] and [D 

• The transverse momenta of the two muons (shghtly smaller for bbnj than for the 
signal, and much larger for the other backgrounds). 

• The distance between them and the closest jet, which is a good discriminator 
against ttnj but not against bbnj. 

• The rapidity difference between the muons and the W* from decay, which is 
reconstructed from the two jets with highest pt- 

• The transverse momenta of the two jets with largest pt- Again, these two vari- 
ables are excellent discriminators against high-p^- backgrounds like ttnj and di- 
boson production, but not very useful for bbnj. 

• The missing transverse momentum. 

• The invariant mass of the two muons and two jets with highest px, ^fifijj- For 
the signal, these four particles result from the decay of an on-shell W boson, 
so the distribution is very peaked around 100 GeV (the position of the peak is 
displaced as a consequence of pile-up, which generates jets with larger pt than 
the ones from the signal itself). Unfortunately, for bbnj the distribution is quite 
similar. 

• The number of b tags and the jet multiplicity. 

• The azimuthal angle (in transverse plane) between the two muons, 0^^. For bb 
without additional jets this angle is often close to 180°, but for bbj and higher 
order processes (which are also huge) this no longer holds. 

The resulting log-likelihood function is presented in Fig. [HI As it can be easily noticed 
with a quick look at the variables presented, the kinematics of bbnj is very similar to 
the signal and so this background is very difficult to eliminate. In particular, for larger 
mjsf requiring large transverse momentum for the leptons drastically reduces bbnj (as 
seen in the previous subsection), but for m^r < it reduces significantly the signal 
as well. As selection cut we require log^g-^sZ-^s — 2.2 for the three final states, 
which practically eliminates all backgrounds except bbnj. The number of remaining 
background events is given in the right part of Table H] (numbers of background events at 
pre-selection equal those in Table [T], and are quoted on the left for better comparison). 
Requiring larger Ls/Lb for the e^e^jj and fJ^^e^jj channels does not improve the 
results, because it decreases the signals too much. The resulting statistical significance 
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Figure 7: Normalised distributions of several discriminating variables for the m^r 
GeV signal and the three background classes (see the text). 
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GeV signal and the three background classes (see the text). The last plot corresponds 
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for the heavy neutrino signals are collected in Table [5l assuming a 20% systematic 
uncertainty in the backgrounds. From these significances, the following limits can be 
extracted: 

(a) A 60 GeV neutrino coupling only to the muon can be discovered for mixings 
|V;.iv| > 0.0070, and bounds \Vf,N? < 1.65(1.95) x 10"^ can be set at 90% (95%) 
CL if a background excess is not observed. These figures are ~ 35 times worse 
than in previous parton-level estimates which overlooked the main background 
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Pre-selection Selection 







e±e± 






e±e± 




TV (a) 


427.3 








42.1 








N (b) 





174.7 








33.9 





iV(c) 


214.0 


88.5 


290.9 


20.4 


17.1 


39.3 


bbnj 


14800 


52000 


82000 


10.7 


291 


96 


ccnj 


(11) 


300 


200 


(0) 








ttnj 


1162.1 


8133.0 


15625.3 


0.3 


1.3 


1.3 


tj 


60.8 


176.5 


461.5 


0.0 


0.0 


0.1 


Wbbnj 


124.9 


346.7 


927.3 


0.2 


2.4 


1.3 


Wttnj 


75.7 


87.2 


166.9 


0.0 


0.0 


0.0 


Zbbnj 


12.2 


68.9 


117.0 


0.0 


1.4 


0.2 


WWnj 


82.8 


89.0 


174.8 


0.0 


0.0 


0.0 


WZnj 


162.4 


252.0 


409.2 


0.6 


0.4 


0.5 


ZZnj 


3.8 


13.3 


12.9 


0.0 


0.5 


0.1 


WWWnj 


31.9 


30.1 


64.8 


0.9 


0.0 


0.0 



Table 4: Number of i^i^jj events at LHC for 30 fb~^, at the prc-sclcction and selection 
levels. The heavy neutrino signal is evaluated assuming m^r = 60 GeV and coupling 
(a) to the muon, V^n = 0.01; (b) to the electron, V^n = 0.01; (c) to both, V^n = 0.01 
and Vi,N = 0.01. 

/i'*'^^ e^e''^ /x^e^ 
N (a) lO.OcT 
N (b) - 0.54cT 
N (c) 4.9(7 0.28(7 1.75(7 

Table 5: Statistical significance of the heavy neutrino signals in the different channels, 
for a mass uin — 60 GeV and couphng (a) to the muon, — 0.01; (b) to the electron, 
VeN = 0.01; (c) to both, VeN = 0.01 and V^n = 0.01. 

bbnj, but would still improve the direct limit from L3 [32,33] by an order of 
magnitude. 

(b) A 60 GeV neutrino coupling only to the electron can be discovered for mixings 
\VeN\ > 0.030, and bounds l^jvp < 3.1(3.6) x 10"^ can be set at 90% (95%) GL 
if a background excess is not observed. 
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The general limits for a heavy neutrino coupling to the electron and muon are displayed 
in Fig. [HI It is interesting to observe that the direct limit from non-observation of like- 
sign dileptons at LHC will be more restrictive than indirect ones from fi — e LFV 
processes at low energies. 

14 
12 
10 

X 

4 
2 

% 0.01 0.02 0.03 0.04 

Figure 9: Combined limits on VeN and V^n, for V^^r = and rrijy = 60 GeV. The red 
areas represent the 90% CL limits if no signal is observed. The white areas correspond 
to the region where a combined statistical significance of 5a or larger is achieved. The 
indirect limit from /i — e LFV processes is also shown. 




4.3 Opposite-sign dilepton signals 

In final states e^fi'^jj the analysis is similar but the backgrounds are much larger. In 
particular, opposite-sign lepton pairs from bbnj production are much more abundant 
than like-sign pairs. Opposite-sign dileptons are produced from tinj dileptonic decays 
and W~^W~nj production (which is larger than W^W^nj). We assume a heavy Dirac 
neutrino with a mass of 60 GeV and VeN = 0.01, V^tv = 0.01. A Majorana neutrino 
gives this signal too, but with half the cross section for the same couplings. We use 
the same pre-selection cuts as in the like-sign dilepton analysis but requiring instead 
opposite charge for the leptons. The number of signal and background events at pre- 
selection is collected in the left column of Table El At pre-selection the bbnj, tinj and 
WWnj backgrounds are 7, 15 and 70 times larger, respectively, than the corresponding 
ones for ix^e^ in Table |H 

The kinematical variables useful for discriminating the signal against the back- 
grounds are the same as for a 60 GeV heavy Majorana neutrino in the like-sign dilepton 
channels. However, in this case the distributions for some backgrounds, namely tinj 
and WWnj, are different. We have grouped backgrounds in three classes: bbnj, tinj. 
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Figure 10: Normalised distributions of several discriminating variables for a 60 GeV 
Dirac neutrino and the three background classes (see the text). 
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Figure 11: Normalised distributions of several discriminating variables for a 60 GeV 
Dirac neutrino and the three background classes (see the text). The last plot corre- 
sponds to the log-likelihood function. 
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Pre-selection 


Selection 


N (e,/i) 


593.5 


14.7 


bbnj 


602000 


73 


ccnj 


5750 





ttnj 


233135.1 


0.3 


tj 


1003.8 


0.0 


Wbbnj 


927.6 


0.0 


Wtinj 


197.0 


0.0 


Zbbnj 


180.8 


0.0 


WWnj 


12016.5 


0.7 


WZnj 


412.1 


0.0 


ZZnj 


14.2 


0.0 


WWWnj 


131.4 


0.0 



Table 6: Number of fi^e^jj events at LHC for 30 fb~^, at the pre-selection and selection 
levels. The heavy neutrino signal is evaluated assuming rriN = QO GeV and coupling 
to electron and muon VeN = 0.01, V^at = 0.01. 

and the other backgrounds (dominated by WWnj). The distributions for the relevant 
kinematical variables and the log-likelihood function are collected in Figs. [10 and [Til 
For event selection we require logj^Q-L^/Lg > 2.9, yielding the number of events in the 
right columns of Table [HI The significance of the heavy Dirac neutrino signal is only 
0.86(7. The combined limits on VeN and V^^ are presented in Fig. [T21 The shape of the 
regions is drastically different from Figs. [6] and [9] because the sensitivity in the e~^e~jj 
and fJ'~^fi~jj channels is negligible, and only when couples sizeably to both electron 
and muon the heavy neutrino signal is statistically significant in the fJ^^e^jj channel. 
The direct limit from non-observation of a fJ^^e^jj excess has a similar shape as the 
indirect limit but it is less restrictive in all cases. 

5 Estimates for Tevatron 

The observability of heavy neutrino signals in like-sign dilepton channels at Tevatron 
seems to be dominated by the size of the signal itself. In contrast with LHC, back- 
grounds are much smaller. For example, the WZjj and WWjj backgrounds have 
cross sections of 0.1 and 0.09 fb, respectively, with the cuts in Eq. (Q. Then, it seems 
reasonable to estimate the total background for 1 fb~^ (including bb) as one event. This 
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Figure 12: Combined limits on VeN and V^a? for a 60 GeV Dirac neutrino. The red 
areas represent the 90% CL hmits if no signal is observed. The white areas correspond 
to the region where a combined statistical significance of 5a or larger is achieved. The 
indirect limit from n — e LFV processes is also shown. 

rough estimation is in agreement with the detailed calculation in Ref . [34] , in which bb 
is estimated from real data. Therefore, if signal events have not been observed with the 
already collected luminosity, upper limits of 3.36 and 4.14 events [35] can be set on the 
signal, at 90% and 95% CL, respectively. From Fig. O and for a fixed mass uin = 60 
GeV, this implies upper bounds \V^n\'^ < 1.3 x 10"^ (90% CL) IV^n]"^ < 1.6 x 10"^ 
(95% CL). This would slightly improve the limits from L3 [32,33]. Of course, a de- 
tailed simulation with the already collected data is necessary to make any claim, and 
the limits eventually obtained will depend on the actual number of observed like-sign 
dilepton events. 

Note also that, given the cross sections in Fig. [21 for a luminosity of 1 fb~^ and 
Vfj,N = 0.098 the heavy neutrino signals only exceed a handful of events for heavy 
neutrino masses uin < Mw, and thus the Tevatron sensitivity (when acceptance and 
efficiency are taken into account) is limited to this mass range. This also holds for 
a heavy neutrino mixing with the tau lepton, for which the N production can be 
larger but r decay branching fractions must be also included in the final cross section. 
Then, if the small excess found by CDF [28] is confirmed, its explanation through 
heavy neutrinos requires additional interactions, for example mediated by a new Z' 
boson [11]. 



24 



6 Conclusions 



Large hadron colliders are not in principle the best place to search for new heavy neutral 
leptons. However, Tevatron is performing quite well and LHC will start operating 
soon, so one must wonder if the large electroweak rates available at large hadron 
colliders allow to discover new heavy neutrinos, given the present constraints on them, 
or improve these constraints. This is indeed the case in models with extra interactions 
[7,8,10]. In this work we have, however, assumed that no other interactions exist and 
that heavy neutrinos couple to the SM particles through its small mixing with the 
known leptons. 

Heavy Dirac or Majorana neutrinos with a significant coupling to the electron can 
be best produced and seen at e^e^ colliders in e^e^ — > Nv iW which has a 
large cross section and whose backgrounds have a moderate size [15, 17,25,36]. On 
the contrary, a Majorana N mainly coupling to the muon is easier to discover at a 
hadronic machine like LHC, in the process qq' — > ^'^N with subsequent decay 

— > ^'^W — > fi'^qq' (plus the charge conjugate). However, even this LNV final state 
is not easy to deal with. SM backgrounds are large and require a careful analysis, 
especially those involving b quarks, for example bbnj and tinj which are the largest 
ones. 

For the simulation of the i^i^jj signals from heavy neutrinos we have implemented 
heavy neutrino production in the ALPGEN framework. In the fi^fi^jj channel we 
have shown, using a fast detector simulation, that a hevy neutrino with a mixing 
Vfj,N = 0.098 can be discovered with a 5a significance up to masses itln = 200 GeV. In 
the region rriN < My/ we find that a 60 GeV neutrino can be discovered for mixings 
|V^7v| > 0.0070; upper limits \V^n? < 1.65(1.95) x 10^^ can be set at 90% (95%) 
CL if a i-i'^i-i'^jj excess is not observed. These figures are in sharp constrast with 
previous estimates, and correspond to the increase in the background estimation of 
about two orders of magnitude (three for mjv < Mw)- In particular, special care 
has to be taken with bb plus jets. The probability of a bb pair to give two like-sign 
isolated muons is tiny, but on the other hand the bb cross section ~ 1 //b is huge. A 
reliable background calculation requires solving this • oo indetermination, what is a 
computationally very demanding task in which some apparently reasonable simplifying 
assumptions, like requiring high transverse momenta of b quarks at generator level, can 
result in an underestimation by a factor of 30. The bbnj background has been found 
to be negligible for larger mjv values but dominant for mjv < M^y (after cuts, 5 times 
larger than the sum of the other backgrounds). This behaviour is due to the very 
different signal kinematics in these two cases. For ttin < Mw the charged leptons 
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are produced with very small transverse momentum, therefore a cut on this variable, 
which could be efficiently used to remove bbnj, cannot be applied. On the other hand, 
requiring e.g. that one charged lepton has px > 30 GeV and the other one Pt > 20 
GeV hardly affects the signal for mjv = 150 GeV, while it practically eliminates bbnj. 

For the other like-sign dilepton channels, e^e^jj and //^e^jj, the prospects are 
worse because backgrounds are much larger. We have found that a heavy neutrino 
with VeN = 0.073 can be discovered with 5a up to masses rriN = 145 GcV. In the 
region mjv < M^, a heavy neutrino with m^r = 60 GeV can be discovered for mixings 
IKjvI > 0.030; upper limits IKatP < 3.1(3.6) x 10"^ can be set at 90% (95%) CL if a 
background excess is not observed. The latter limits are of the same magnitude but 
worse than those from L3. Besides, couplings of this size would be in conflict with the 
non-observation of neutrinoless double beta decay, requiring cancellations with other 
new physics contributions. Finally, for a heavy neutrino with mjy = 60 GeV and 
coupling to both electron and muon we have found that direct limits on Vg^ and V^jy 
will improve the indirect ones from fi — e LFV processes. For completeness we have also 
examined the LHC sensitivity for a Dirac neutrino coupling to the electron and muon, 
in n^e^jj final states. The sensitivity is much worse, as expected from the larger LNC 
backgrounds involving opposite-sign dilcptons, and the direct limits obtained are worse 
than the present indirect ones. Hence, LHC is not expected to provide any useful direct 
limit on heavy Dirac neutrinos, for which all final states conserve lepton number. By 
the same token, other decay channels such as N ^ Zi/, N — > Hi/ and production 
processes as pp ^ Z ^ Nv, have too large backgrounds as well. 

In the detailed analyses presented for itln = 150 GeV and mjv = 60 GeV we have 
shown that background suppression {ttnj and diboson production in the former case, 
bbnj in the latter) is not efficient with simple kincmatical cuts, and requires more 
sophisticated methods, like the likelihood analysis applied here, or neural networks. 
The analysis could be further improved when one includes other variables not accessible 
at the level of fast simulation. For example, a bb pair giving two like-sign isolated muons 
is most often caused by the oscillation of one of the mesons before decay. This 
should appear as a secondary vertex, which could be identified in the detector. On the 
other hand, the possibility of lepton charge misidentification should be addressed. Full 
simulation of bbnj for the LHC luminosity is beyond reach of present and foreseable 
computers, so this background will have to be estimated from data. In any case, we 
stress that bbnj, as well as ttnj, must always be considered as a potentially dangerous 
source of two like-sign dileptons. And, if a moderate background excess is observed 
at low transverse momenta, a precise evaluation of the bbnj background is compulsory 
before drawing any conclusion. 
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It is finally worth noting that heavy neutrino decays, as for any other vector-like 
fermion, are a source of Higgs bosons [37]. Nevertheless, in contrast with the quark 
sector [38] Higgs boson production from N decays is rather small, and only a handful 
of fiN fiuH jivbh events are expected to be found at LHC. Besides, we also 
point out that large effects due to heavy neutrinos and, more generally, other neutrino 
physics beyond the SM might be observed at large hadron colliders. However, in 
all cases they require new interactions and often model dependent constraints. This 
means further assumptions, and in this situation the main novel ingredient is not 
only the heavy neutrino. In this category there are many interesting scenarios, also 
including supersymmetry (see for an example Refs. [39,40]). Then, compared to these 
new physics models the limits established in this work are modest. For example, if 
the heavy neutrino has an interaction with a typical gauge strength, as in left-right 
models with a new heavy Wr, the LHC reach for tjin increases up to approximately 
2 TeV [7,8]. In the case of a new leptophobic Z'^ boson, the mjsi reach in N pair 
production pp ^ Z'^ ^ NN is increased up to 800 GeV [11]. 
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A Evaluation of the bb background 

bb production, which has a huge cross section of order 1 /ib at LHC, is the largest SM 
source of like-sign dileptons. Charged leptons are produced in the decays b ci'iy, 
b —>■ ci'^u, and like-sign lepton pairs can arise when one of the b quarks yields a 
5° meson which oscillates before decay. Additionally, like-sign charged leptons can be 
produced from the subsequent decay of the charm (anti)quark, e.g. c si'^u. We have 
investigated the relative contribution of the two sources by simulating with Pythia a bb 
sample of 25 million events with and without mixing. The number of dielectron and 
dimuon events (requiring isolation and transverse momentum greater than 10 GeV) is 
gathered in Table [71 A quick look at these numbers reveals that about 80% of like-sign 
dileptons result from oscillation. 
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B mixing No B mixing 





55 


12 




456 


109 


e±e± 


1242 


334 




309 


335 




1357 


1643 


e±e± 


3755 


4671 



Table 7: Number of dilepton events obtained from a sample of 25 million hh events, 
when B^ oscillation is allowed in Pythia (first column) or not (second column). 

It must be emphasised that the theoretical evaluation of the hh contribution to the 
like-sign dilepton SM background involves several uncertainties. The most obvious one 
affects the total hh cross section, which depends to a large extent on the generation 
cuts placed on h transverse momenta. A second one involves h quark fragmentation. 
We have used the Peterson parameterisation with = 0.0035 [41]. With the default 
Pythia setting = 0.005 the number of (isolated) dileptons obtained is a factor ~ 0.77 
smaller. But perhaps the largest uncertainty comes from the fact that our analysis relies 
on a fast simulation of the detector, which may be inadequate when studying delicate 
issues like lepton isolation. At any rate, a full simulation of a large bb sample is out of 
present reach and this background will have to be measured using real data. 

Apart from these theoretical uncertainties there is another one due to the limited 
statistics of the simulated samples. The hh cross section is 1.4 fih when both h quarks are 
required to have > 20 GeV at the generator level. Fast simulation of 30 fh^^ would 
take about 15000 days in a modern single-processor system, making this computation 
difficult even in multi-processor grids. (Full simulation would take about 10^ years and, 
as emphasised above, in the real experiment this background must be estimated from 
data, as it has been done by DO [34].) Therefore, for our evaluations we have simulated 
samples of approximately 100, 35, 15 and 5 million events for bb, bbj, bb2j and bbSj, 
respectively, corresponding to a luminosity L — 0.075 fb~^ and the cross sections given 
by ALPGEN. The size of the samples is reduced when performing the MLM matching, 
which has efficiencies of 90.7%, 41.8%, 18.7% and 12.7%, respectively. The number of 
events at pre-selection is calculated by rescahng the number of events to 30 fb~^. For 
example, 

N(i,^l,^; pre, 30) ~ A^(/x±//±; pre, L)/i , (7) 

with fi — 400. This rescahng introduces a large statistical uncertainty and, moreover, 
the estimation of the number of events after selection cuts cannot be done in this 
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iV(/i±/i±; sel, 30) ~ iV(££'; sel, L) 



-h 



(8) 



way, since no fi^fi^ events survive the cuts applied. Instead, we make the reasonable 
assumption that selection cuts, which are based on kinematical variables, have the 
same effect on all if events, where i, i' = e, /i, not necessarily of the same charge. 
Then, for bbnj backgrounds the number of fj^^fi^ events after selection cuts for 30 fb~^ 
can be estimated from the samples with a smaller luminosity L as 

A^(/i^/i^; pre, L) 
N{ee';pTe,L) 

Since the total number of ii' events is about 200 times larger than the number of /x^/i^ 
events, the term in brackets in Eq. ([8]) is of order two, and thus the simulated samples 
provide a statistically more precise estimate of the results for yU^/x^ final states. We 
have explicitly checked whether the relevant kinematical distributions are similar or 
not for several final states. In particular, differences between electrons and muons 
might be expected due to the different energy resolution and isolation criteria. The 
most crucial variables for background suppresion are the transverse momenta of the 
two leptons. They are presented in Fig. [121 together with a "signal" sample included 
for comparison. The heavy neutrino sample corresponds to more than 40000 events, 
while the dilepton samples from bbnj contain 37 fj^^fi^, 382 fj^~^fi~ , 1497 e^e^ and 4676 
e'^e~ events, respectively. The fi'^fi~ , e^e^ and e"^e" samples have remarkably similar 
distributions, while fi^fi^ events apparently concentrate at lower transverse momenta. 
This seems to be only a statistical effect, given the smallness of the sample (only 37 
events). This belief is strengthened if one realises that e^e^ and e~^e~ events have 
the same distributions (what suggests charge independence) and the same happens 
for e~^e~ and fi'^fi~ (suggesting flavour independence). Two further variables which 
might exhibit differences are the distance between the leptons and the closest jet, also 
shown in Fig. [131 In this case there seem to be small differences between the samples. 
However, these two variables are not determinant in suppressing the background, as 
it can be observed by comparing with the signal distribution, and any eventual 
difference in kinematics will have little effect on our calculations. 

Finally, it is worth remarking here that raising the threshold at generator level, 
e.g. to 50 GeV, leads to a dramatic reduction of the bbnj cross sections, making the 
simulation more manageable. However, this also results in a gross underestimation of 
the bbnj background. We have checked this by simulating two samples of 25 million bb 
events with > 20 and p^ > 50 GeV, respectively. For pre-selection we just require 
two isolated muons of either charge with pf^ > 10 GeV. For the p^ > 20 sample we 
obtain 364 bb —* fifi events, while for the > 50 sample we only obtain 287 events. 
Given the difference in cross sections (1430 nb for > 20 GeV and 58.8 nb for p^^ > 50 
GeV) , raising p^ to 50 GeV at event generation would underestimate this backgrounds 
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Figure 13: Distribution of several kinematical variables for /i^/x^ events from heavy 
neutrino production and dilepton events from bb production (see the text). 



by a factor of 30. For bj we have checked that raising p\, at event generation to 50 
GeV reduces the number of dimuon events by a factor of 35. This seems to be related 
to the fact that b quarks with larger transverse momentum give more collimated decay 
products, and thus the muons are less isolated. On the other hand, b quarks with 
too low transverse momentum cannot produce muons with pi^ > 10 GeV. For this 
reason, we expect that our evaluation of bbnj provides a good estimate of the actual 
background to be found at LHC. 

The evaluation of from ccnj proceeds in the same way. However, the number of 
dilepton events is much smaller and no /i^/x^ events appear in the samples simulated 
(containing about 145 million events after MLM matching). In this case the number 
of /i^/i^j j events from ccnj production is estimated by comparing with e^e^nj events 
and assuming the same ratio of events N{fi^fi^jj)/N{e^e^jj) as for bbnj production. 
The result is shown between parentheses in the tables. 
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B Heavy neutrino mass reconstruction 



For heavy neutrinos heavier than the W boson the decay — > fiW fiqq' takes 
place with W on shell; thus, the invariant mass of the two quarks is up to finite 
width effects. In simulated signal events, however, several extra jets often appear due 
to pile- up and initial/final state radiation, and it is not straightforward to identify the 
two ones originating from the W decay. We have tested two procedures: 

(i) To take, naively, the two jets with highest transverse momentum. This method 
will be denoted as 'Rl'. 

(ii) To try all possible pairings among the jets, choosing the pair giving an invariant 
mass closest to M^/. This method is denoted as 'R2'. 

The difference between the two choices is illustrated in Fig. [13] (left) for the case of the 
heavy neutrino signal. The method Rl yields a moderate peak for the W reconstructed 
mass M^*^. When included in likelihood function (see Fig.H]), this variable improves the 
signal significance about 2%. (No improvement is found when performing a kinematical 
cut on in addition to the cut on likelihood.) The second method R2 gives a 

considerably more peaked distribution for the signal, at the expense of strongly biasing 
the background, as it is shown in Fig. [H] (right). Thus, with the second method the 
W invariant mass is not a useful variable for discriminating the signal against the 
background. 
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Figure 14: Left: reconstructed W mass for the heavy neutrino signal, using the two 
methods (Rl and R2) explained in the text. Right: reconstructed W mass for the 
signal and SM backgrounds using the method R2. 



The heavy neutrino mass is obtained as the invariant mass of the jet pair selected to 
reconstruct the W, plus one of the two muons. In order to improve the reconstruction. 
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the two jet momenta are rescaled so that their invariant mass coincides with M^. For 
both W reconstruction methods the results are very similar, as it can be observed in 
Fig. [151 The invariant mass of the W and the muon with smaller transverse momentum 
'nT'Wfi2 is niore concentrated around the true rriN, and is taken as the heavy neutrino 
reconstructed mass in our analysis. In case of discovery, this distribution might be 
used to determine rriN- 



0.15 




0.15 



200 



iV(150)Rl 
A'(150)R2 




Figure 15: Invariant mass of the reconstructed W boson (rescaled) and the muon 
with highest (/xi) and lowest {fi2) transverse momentum, for the two W reconstruction 
choices explained in the text. 



For rriN < Mw it is very difficult to identify the two jets coming from — > W*jj,, 
which have low transverse momenta, due to the appearance of extra jets from pile-up. 
This fact is clearly seen examining the invariant mass distribution of the two jets with 
highest pt and either of the two muons, in Fig. [161 In both cases the distribution peaks 
well above mjv = 60 GeV, indicating that one or the two jets do not really originate 
from the heavy neutrino decay. We have not found any improvement of the signal 
significance considering these variables in the likelihood analysis. 

In case of discovery, one possibility for the mass determination could be to 
consider the invariant mass of the two muons, which we present in Fig. [T7] (left) for 
heavy neutrino masses of 50, 60 and 70 GeV. This distribution seems to peak around 
mAr/2. Other possibility is to exploit the fact that, since the on-shell decay W —>■ fiN 
is two-body, the energy of this muon in the centre of mass (CM) system, -E^^, is fixed 
by rriTv- Thus, we may determine the heavy neutrino mass as 

m^M = - 2M^ECM . (9) 

The mjv reconstruction from the muon energy in the CM frame (defined as the rest 
frame of the two muons and two jets with largest pt) is shown in Fig. [TT] (right). For 
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— Af(60)n, 




Figure 16: Invariant mass of the two jets with highest transverse momentum and each 
of the two muons, for a heavy neutrino signal with niN = GO GeV. 



each event two values of m'^ are calculated, corresponding to the two possible muon 
choices, and both are plotted. Imaginary values are discarded. These procedures for 
mjv determination will be subject to possibly large systematic uncertainties, but their 
evaluation is beyond the scope of this work. (For example, the reconstruction from the 
muon energy in the CM frame is expected to have a systematic uncertainty from pile- 
up, which could be decreased using the muon energy in the laboratory frame, but at 
the expense of losing sensitivity to niN-) If heavy neutrinos were discovered, interesting 
information about CP violation, relevant for leptogenesis, could also be inferred [42]. 
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Figure 17: Left: invariant mass of the two muons, for three heavy neutrino masses. 
Right: heavy neutrino mass reconstructed from the muon energy in CM frame. 
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